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Abstract 

The Stieltjes constants jk( a ) appear as the coefficients in the regular part 
of the Laurent expansion of the Hurwitz zeta function £(s, a) about s = 1. We 
present series representations of these constants of interest to theoretical and 
computational analytic number theory. A particular result gives an addition 
formula for the Stieltjes constants. As a byproduct, expressions for deriva- 
tives of all orders of the Stieltjes coefficients are given. Many other results are 
obtained, including instances of an exponentially fast converging series repre- 
sentation for 7^ = 7fc(l). Some extensions are briefly described, as well as the 
relevance to expansions of Dirichlet L functions. 
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Introduction and statement of results 

The Stieltjes (or generalized Euler) constants 7* (a) appear as expansion coeffi- 
cients in the Laurent series for the Hurwitz zeta function ((s, a) about its simple pole 
at s = 1 [3 [71 [TH [231 12Z] - These constants are important in analytic number theory 
and elsewhere, where they appear in various estimations and as a result of asymptotic 
analyses. They are also of much use in developing a binomial sum S^(n) introduced 
by the author in the study of a critical subsum in application of the Li criterion for 
the Riemann hypothesis [9|. The constants 7fc(l) are important in relation to the 
derivatives of the Riemann xi function = tt~ s ^ 2 T(s/2 + — 1)C( S ) at s = 1, 
where F is the Gamma function and ((s) is the Riemann zeta function [T^, [T9 | [20 ], [28] . 
Hence some of the relevance to the Li criterion. 

Despite the fact that many relations are known for the Stieltjes constants, there 
are many open questions, including those concerned with their arithmetic nature. 
Even good estimations of their magnitudes is still lacking. In addition, formulas for 
obtaining the Stieltjes constants to arbitrary precision remain of interest. 

On the subject of the magnitudes |7fe(a)|, known estimates pfj [36] are highly 
conservative, and there is much room for improvement. We present one approach 
for such estimation. Our result is probably less important by itself, than for the 
possibilities for extension and improvement that it suggests. 

In this paper, we present various series representations of the Stieltjes coefficients. 
Most of these are fast converging, so have application to high precision computation. 
In particular, we make use of the Stirling numbers of the first kind s(k,j) and their 
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properties. Therefore, our treatment reflects a fusion of analytic number theory and 
enumerative combinatorics. Moreover, one of our Propositions (9) gives instances of 
an exponentially fast converging series representation for 7*.. Such very fast converging 
series are relatively rare, and this may be the first result of its kind for the Stieltjes 
constants. 

The Hurwitz zeta function, initially defined by a) = J2^=o( n + a )~ s f° r Re s > 

1, has an analytic continuation to the whole complex plane [H, [321 ED]- In the case 

of a = 1, C( s ) a ) reduces to the Riemann zeta function ((s). In this instance, by 

convention, the Stieltjes constants 7^(1) are simply denoted 7^ [5l [T8l l23~t |25 | [271 [36] . 

We recall that 70(0) = —if) (a), where ip = T'/T is the digamma function. We also 

recall that jk(o> + 1) = 7^(0) — (ln fc a)/a, and more generally that for n > 1 an integer 

, x ^ln fc (a + j) , . 

lk{a + n)= lk a - £ \ J> , 1.1 

as follows from the functional equation £(s, a + n) = £(s, a) — X^=d ( a + j) s - 

Unless specified otherwise below, letters j, k, £, m, and n denote nonnegative inte- 
gers. We point out the common alternative notations and ^ for Stirling num- 
bers of the first kind s(n,m) [TJ [131 H7J) with the relation '" = (— l) n+m s(n, m) 



n 
m 

n 
m 

holding. We also make use of P%(t) = B\(t — [t]) = t — [t] — 1/2, the first periodized 



Bernoulli polynomial (e.g., [19, 32J). Obviously we have |-Pi(t)| < 1/2. 

Proposition 1. (Addition formula for the Stieltjes constants) Let Re a > and 
\a\ > \b\. Then (i) 

le (a + b)= 7i (o) + (-l) e £ JJ—^ EC" 1 )* (I) k + mC^U, a). (1.2) 



(ii) We have 

7 ;(a) = (-l) m [C W (2,a)+£C ( ^ 1) (2,a)], (1.3) 
7?(a) = (-l)'[2C w (3,a) + 3«C ( '- 1) (3,a)+^-l)C ( '- 2) (3,a)], (1.4) 

and 

7 f (a) = (-l) m [6C W (4, a)+ll£C ( '' 1) (4, a)+6£(£-l)C ( ^ 2) (4, a)+£(£-l)(£-2)C ( '~ 3) (4, a)]. 

(1.5) 

Corollary 1. 

7 ;(l)=C(2)[7 + ln(27r) + 12C / (-l)], (1.6) 

where £(2) = n 2 /6 and 7 = — tp(l) is Euler's constant. 
Corollary 2. For j > 1 we have 



It 



(?) 



'I 



(a) = (-lfE(-l)^ [ k )8{j + l,k+ l)(V- k \j + 1, a). (1.7) 



fc=0 



Proposition 2. (i) For Re a > 1, 

7 „(a) = -A— ^ - n! £ -— - £ M"«(* + 1, n - ro + l)C (m) (* + M)> 

n + 1 fc=i r + i J ! m =o m - 

(1.8) 

(ii) For Re a > 1/2, 



ln n+1 (a- 1/2) ~ 1 " (-l) m 



7»(«) ^-!E^E^#H- + l)('l2 W)fl ) ) 

(1.9) 

(iii) For Re a > 1/2, 

7 „ (a) = J.^V - 1/2) + „, (=g. (a+ l, n- m+ l) C '-'(2. + l, „) 



~ 2n] g 16*(2fc+l)! ^ hr s{4k + 1 > n ~ m + 1 )C M (4A: + 1, a). (1.10) 
Proposition 3. (Asymptotic relation) Let Bj denote the Bernoulli numbers. For 
Re a > and I — > oo we have 

7<(a) . i. ^ „ _ *g± _ £ _^„- £ Q Hs(2m , t + i) ln « „. (1 .n) 

Proposition 4. For any integer N > and Re a > we have 



£!^j ! ^g^^ ( ^ riW)[( _ 1) ^ (ri . 



v i^-fc 



ln *(n + o) 



;i.i2) 



Proposition 5. Let a > 0. Then there exists a* > 1 such that 71(a) is mono- 
tonically increasing for a < a* and 71(a) is monotonically decreasing for a > a*. The 
approximate numerical value of a\ is a* ~ 1.39112. 

Proposition 6. Let Re z > 0. Then we have 

1 00 z k r 1 (— IV 

^m! 11 *^' (L13) 

Proposition 7. Put as in [251 [36] , C n {a) = j n (a) - \ ln n a for < a < 1. Then 

we have 

|Cn(«)l<77^ ; (L14) 

Our methods extend to many other analytic functions. Additionally, although we 
present explicit results for expansions at s = 1, this is not a restriction. For those 
functions possessing functional equations, we also gain expansions typically at s — 0. 



More generally, expansions about arbitrary points in the complex plane are usually 
possible. 

Another instance for which our methods apply is for the Lerch zeta function 
$(z, s, a). As an example, we consider the Lipshitz-Lerch transcendent 

oo 2ninx 

L(x,s,a) = y- = <S>(e 2wix ,s,a), (1.15) 

1 J ^ (n + a)* 

for complex a different from a negative integer. We take in (1.15) x real and non- 
integral, so that convergence obtains for Re s > 0. Else, for x an integer in (1.15), 
we reduce to the Hurwitz zeta function. The functions $ and L possess integral 
representations and functional equations. 

A case of particular interest for (1.15) is when x = 1/2. Then we obtain the 
alternating Hurwitz zeta function, 

Therefore, expansion at s = 1 yields expressions for differences of Stieltjes constants 
7fc(a/2) — 7fc[(a + l)/2]. More generally, for x Z, L is nonsingular at s = 1 and we 
may write 

oo / -| \n 

L(x,s,a) = ]T ^-J-£ n ( x ,a)(s - l) n . (1.17) 

n=0 n - 

We have 

Proposition 8. Let a > and |£| < a. Then we have the addition formula 

oo tk-\ oo / \ / . \3 

£ n (x, a+0 = tn{x, a)+E ^~ iy [jj ( n -^< k i {ds) HX ' S + k ~ 1j a)ls=1 

(1.18) 



As a result, we obtain for derivatives with respect to a of £ n , 
Corollary 3. 

a) = fX-iy Q (n - j>(* + 1, n - j + 1) (±J L(x, s + k, a)\ s=1 . (1.19) 

Let r(a, #) = J^ 00 t a_1 e _i <it be the incomplete Gamma function, p F g the generalized 
hypergeometric function, Ei(z) = — J*^ e _< (c?t/t) the exponential integral, and erf the 
error function [H [16]. As a foretaste of a family of other results we offer 
Proposition 9. We have (i) 

OO 1 OO I 

" n=l n=l 

and (ii) 



X) -[1 - erf(V5m)] - £ Ei(-7m 2 ) - 1 + - 1h(4tt) 



;i.20) 



1 OO 1 

-V- 

2^" 



7x = — + -tb (-} + -ib 2 - 1 + - In vr - - In 2 vr 
' 16 T \2J 8^ \2J 2 8 

4n 2 F 2 (-, -;-, -; -n 2 7r) + ^ - 2 In n - erf (ny^) hi 7r 



00 r 7 2 , vr 2 1_ 2 _ 



+ X | X + 24 " i n2?F 3Fs(1 ' 1; 2 ' 2 ' 2; ""^ + I ln(n27r) [27 + ln( ™ 2)] + ^ ln ^(-^O 

(1.21) 

where V(V 2 ) = -7- 2 In 2. 

We introduce the polylogarithm function, initially defined by Li s (z) = Y^kLi f° r 
|z| < 1 and Re s > 1, and analytically continued thereafter. We illustrate a method 
that more generally leads to expressions for the sums 7&(a) + 7&(1 — a). We have 
Proposition 10. Let < a < 1. We have (i) 



1\ 9 
In ir+ib ( — J — 7r cot Tca — 2'ip(a) = 7+ln7r + 2 — 



[Li s (e 27rM ) + Li s (e- 2 ™)], (1.22) 



s=0 



and (ii) 



7i( a )+7i(!-«) = tt^ 



vr 2 ln 2 vr lnvr , /1\ 1 ,„ /1\ 1 



1 



12 4 2 

d 



2) 4 r \2J 2 



+ - 



-^■(7+ln7r) 2 -(7+ln7T. 
4 os 



[Li s (e 27rm )+Li s (e- 27rw )] + 



s=0 



In 7r — ip 



(a)H-V'(l-a)] 



[Li s (e 2 ™)+Li s (e- 



-27ria\ 



;i.23) 



Lastly, we present expressions with Stirling numbers for rapidly converging ap- 
proximations to 7 = 7o- Historically this subject [29J has been important especially 
due to P. Appell's use of them in his attempted proof of the irrationality of Euler's 
constant [3]. In addition, we present an exact representation for the difference be- 
tween the Stieltjes constants and a finite sum. Near the end of the paper, we further 
discuss these subjects. We have 
Proposition 11. Define polynomials for n > 1 



1 [V 

P n +i{y) = -r / x(l - x)(2 - x) ■ ■ ■ (n - 1 - x)dx, ( 
n\ Jo 



1.24) 



and their values p n +i = Pn+iO-)- Put r n^ — lk~ where 



D (k) = y 



In m 



ln fc+1 (n + 1). 



m=l 



m k + 1 



1.25) 



Then we have (i) 



Pn+M 



' )U E(-l) fc ^7^ l) m-l)y+y+m-y) k ~\y-l)+ll (1.26) 



k=l 



k(k + l) 



Pn+M 



\n+l n 



in j k+1 



(-1) 



n+1 



n-l 



s(n, fc) 



n + 1 



(1.27) 



where Sjk is the Kronecker symbol, (iii) the special case 



fe=i 



and (iv) 



where 



m=n+l 



(1.28) 



(1.29a) 



OO J 



fc-1 



A:! 



(-l) n ln fc m + — J] 71 — -rr;s(n + 1, j + 2) In* - ' -1 m 



(1.296) 



Proof of Propositions 

Proposition 1. We let (z)k = T(z + k)/Y(z) be the Pochhammer symbol. We make 
use of the following. 
Lemma 1. We have 



'cT 

ds 



:-i) j+e e\ s (j + i,£ + i). 



(2.1) 



s=l 



We have the standard expansion 

(s)j = s(s + 1) • • • (a + j - 1) = 



(2.2) 



k=0 



giving 



d_ 

ds 



= E(- 1 ) J ' +fes (i' *)*(* -!)•••(*;-€ + l)s fc , (2.3) 
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so that 



ds 



s=l 



k=l 
3 



k=e 

3 



= (-i)'*!E(-i)*«0\*) 

k=i 

It is known that [17J (p. 265) 

E(-i)%^)(Tl = 



(2.4) 



(2.5) 



Therefore, Lemma 1 follows. 

In order to obtain Proposition 1 we apply a formula of Wilton [35] for the Hurwitz 
zeta function, 

00 (—1)3 rfs + j) 

C(s,a + 6) = C(s,a) + E ^ C( g + .?>)^ s ^ x > H < H> ^ a > 0, 



fi j« r( s ) 



along with the product rule 



'cT 
ds 



{s)jt{s + j,a) 



s=l 



E 

fc=0 



'cT 
ds 



(2.6) 



(2.7) 



8=1 



The defining Laurent expansion for the Stieltjes constants is 



J_ ~ (-l) fc 7fc (a) 

- 1 



C(a, a) = — - + E V ^,, ,KV ^ ( g - l) fc , s + 1. 



fc=0 



k\ 



(2- 



Therefore, the Wilton formula (2.6) gives 



_iU oo (_\\k oo / \\3 

E Hi + - x ) fc = E 4r 7fe(a)(s " l)k + S ji (s)jC(s + J ' a)6 "' 



00 ^ 1 \ A; 

E 

fc=0 



J=l 
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\b\ < \a\, Re a > 0. (2.9) 
We take £ derivatives of this equation and use (2.7) and Lemma 1, yielding 

= ("1)S(«)+E^E 7 )(-l) k s(j+l,k+l)k\^- k \j+l,a). (2.10) 

Therefore, Proposition (i) follows. 

For part (ii), we form the limit difference quotient 

i t {a) =\\m\[ 1 ,{a + b)- 1 ,{a)l (2.11) 

giving 

7{(a) = (-1)* p(-l) k Q 5(2, k + l)fc!C^(2, a). (2.12) 

The Stirling number truncates the summation, with s(2, fc + 1) = (— l) fc+1 , k — 0, 1, 
and otherwise is for A; > 2. 
We next have 

7 ;'(a) = lim l[ 7 ,(a + 26) - 2 7 ,(a + b) + lt {a)\ 

= (-1) £ [2C W (3, a) + 3£C (£ - 1} (3, a) + €(€ - 1)C ( ^ 2) (3, a)]. (2.13) 

Similarly higher derivatives of 7^ (a) may be determined and part (ii) has been shown. 

From part (ii), we have 7i(l) = C'(2) + C(2), wherein by the functional equation 
of the Riemann zeta function, or otherwise, we have the relations 

C'(2) =C(2)( 7 + ln2-121nA + ln7r) 
11 



= C(2)[ 7 + ln(27r)-l + 12C / (-l)], (2.14) 

where In A = 1/12 — ('(—1) and A is Glaisher's constant. Corollary 1 follows. 

For Corollary 2 we simply shift the summation index j — > j + 1 in the Taylor 
series (1.2) and read off the derivatives. Otherwise, we could make use of the forward 
difference operator &- n A"[/](x) = & _n E£=o ft) (-l) n ~ k f(x + kb). 

Remarks. The auxiliary relation (2.5) can be proved in a number of ways. One 
is with induction by using a recursion relation satisfied by s(j, k) and by the binomial 
coefficient. Other methods include the use of integral representations either for the 
Stirling numbers of the first kind or for the binomial coefficient. For a proof using 
boson operators, see [21] (Identity 1). 

It is readily checked that at £ = Proposition 1 (i) yields the identity 70 (a + b) — 
—ifj(a + b). For we have 



00 

7o(a + 6)=7o + E(- & )'C(j + l,a) 
3=1 




= 70(a) + V(a) - if>(a + b) = -V>(a + 6). (2.15) 

Herein we used a standard integral representation for ((s,a) and for the digamma 
function [I] (p. 259) or [IB] (p. 943). 
Similarly, we may write 

00 

7i (a + b)= 71(a) - ST{-b) j [C{j + !,«) + HjCij + 1, a)], (2.16) 
3=1 

12 



where Hj = Yjk=i x is the usual harmonic number. Again, integral forms of this 
relation may be given. 

Our integral representation [TTJ [Proposition 3(a)], 

/x 1 i a- ln fe+1 a 2^ f°° (y/a-i)\n k (a-iy) , n . 1 . 

lk {a = —ln k a-—— + -Re ) y ' "U y, Re a > 0, 2.17 

2a k+1 a Jo (1 + y 2 /a 2 )(e 27T y — 1) 

could also be used to prove Proposition 1 part (ii). 

The approximate numerical value 7^(1) — 0.707385812532 suggests that 7^(1) = 
C'(2) + C(2) can be written as l/y/2 together with a series of systematic correction 
terms. 

Of course, the Wilton formula has built in the relation ^C( s > a ) — ~ S C( S + 1> a )i 
and more generally that 

^Vc(«,o) = (-l) i (*)iC(* + J,o). (2-18) 

Proposition 2. These results are based upon [8] (Proposition 1, parts (i), (ii), and 
(iv)). We have for each indicated domain of a, 
(n — 1 V~ s 1 00 TV e _i_ M 

c(s ' a) = 3 _i " £ JkTW C(s + ^' a) ' Re a > (2 ' 19) 

C ( s a ) = ( a ~ 1 /2) 1 ' s _J_y T }, S + 2k lx (s + 2k, a), Re a > 1/2, (2.20) 
sv ' ; s-1 r(s) ^4 fe (2A; + l)! sv ' ; ' 7 ' v 7 



and 



is _ 2 (2a-l) 1 - , 1 \l^(-l) k T(s + 2k) 



s — i r (s 



- Y l, , \, C(s + 2k, a) 
2^ Q A k (2k + 1)\ SV ' ; 



- 7 — h—, ttC(s + 4fc, a) 

^ 16 fc (4fc + l)! sv ' 7 
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Re a > 1/2. (2.21) 



Lemma 1 immediately carries over to 



( 



d_ 

ds 



) 



p.i 



8=1 



(-i) pj+e e\s( P j + i,£ + i), P >i. 



(2.22) 



We then expand equations (2.19)-(2.21) in powers of s — 1 using the product rule and 
the Proposition follows. In the case of part (iii), we have first moved the k = term 
on the right side of Eq. (2.21) to the left side and multiplied the resulting equation 
by 2. 

Remarks. Part (ii) of Proposition 2 gives the general n case beyond the low 
order instances given explicitly in terms of generalized harmonic numbers H^> in 
Proposition 7 of Ref. [10J. After all, it is well known that s(n + 1,1) = (— l) n n\, 



s(n + 1,2) = (-l) n+1 n\H n , and s(n + l, 3) = {-l) n n\[Hl- H^}/2, where H n = #«• 



This part of the Proposition has also been obtained by Smith [30] . 

Part (i) of Proposition 2 exhibits slow convergence. In contrast, parts (ii) and (iii) 
are very attractive for computation. 

The Maslanka representation for the Riemann zeta function is written in terms 
of certain Pochhammer polynomials Pk{s) [26J. Therefore, it is also possible to write 
expressions for jj from this representation in terms of sums including Stirling numbers 
of the first kind. 

Likewise, the Stark-Keiper formula for ( may be used to develop expressions for 
the Stieltjes constants. For iV > an integer, this representation reads 



Serviceable series for the Stieltjes constants using the Stirling numbers of the first 




(2.23) 
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kind can also be written using the Taylor-series based expressions 

00 (-a) n 

((s, a) =a~ s + J2 HMsU(s + n), \a\ < 1, (2.24) 

and 

C U a + - ) = £ ^( S ) n (2 s+ " - 1)C(« + n), \a\ < 1/2. (2.25) 

In general, Dirichlet L functions may be written as a combination of Hurwitz zeta 
functions. For instance, for \ a character modulo m and Re s > 1 we have 

^xHgf^gxMC^). (2.26) 

For x a nonprincipal character, convergence obtains herein for Re s > 0. Therefore, 
our results are very pertinent to derivatives and expansions of Dirichlet L series about 
s — 1. Especially for real-Dirichlet-character combinations of low order Stieltjes 
constants, Ref. PU] may be consulted. 

Proposition 3. We have the representation valid for Re s > — (2n — 1), 

_1— s n r) 

53 ( " 

fc=l 



C( S ,a) = a- s + ^r + E(^-i^f 



If we take n — > 00 in this equation we in fact obtain an analytic continuation of the 
Hurwitz zeta function to the whole complex plane. We then develop the result, 

1 „1— s oo r> 

C(s, a) = l -a~ s + ^— + EWi-i^ "^ 1 . ( 2 - 28 ) 

& Si - =2 ]■ 
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as a series in powers of s — 1, where simply the factor 

= aT' (-!)'-* ]n e ~ k a. (2.29) 



, \ e-k 
a \ 



ds 



s=l 

We use the product rule for derivatives of the summation term together with Lemma 
1 and find 

l£ (a) -Ya Wa ~ TTT + ^ 1)1 'f" ' f Q \k) S{h k + m ^ ^ (2 - 30) 

Since i?2n+i — for n > 1, the stated form of the Proposition follows. 

Proposition 4- We employ the representation [2] (p. 270) based upon Euler- 
Maclaurin summation and integration by parts, for N > 0, and Re s > —m, with 
m = 1,2,..., 



v i (N + a y~ s ™ ( s ) r v - i 



n=0 



n + a) s s-1 ^(r + 1)! 



N 



((s + r,a) - J2 



n=0 



(n + a) f 



Is) .i 00 r 1 7 i"i+i 

l£V±i J , (2 31) 

(m + 1)! ^./o (n + a + w) s+m+ v 7 

Taking m — > oo and developing the result in powers of s — 1 using Lemma 1 gives the 
Proposition. 

Remarks. In practice in using Proposition 4, there will be a tradeoff in selecting 
iV and the cutoff or otherwise estimating the remainder neglected in the sum over r. 
We expect the convergence to be poor as Re a — > in this Proposition. We anticipate 
that a suitable procedure for computations is to take values with Re a > 1, and then 
to use relation (1.1) as needed. 
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Similarly, we could employ the representation [16] (p. 1073) for integers N > 
and Re a > 

C( s , a) = V . 1 . + (iV + Q i )1 " - s V ^"^ (ft, Re s > 1, (2.32) 

where the integral may be easily expressed in closed form. Or we could use the 
representation for integers N > and Re a > [20] (p- 16) 

N 1 1 / 1 \ 1 ~ s r°° Pi ft) 

c(^) = Et— u + — 7 r + 9 +a W , r7 r w dt > Res>0 - 

^ (n + a) s s - 1 V 2 / Jn+i/2 (t + a) s+i 

As we easily have that the integral in this equation is bounded by 

1 

N + l/2 + a| s ' 

the integral converges uniformly for s in any compact subset of the half plane Re 
s > (and for arbitrary a). 

Proposition 5. The function 71(a) — > —00 as a — > + and as a — ► 00. Indeed, 
the asymptotic form as a — > 00 is 71(a) ~ — ~ In 2 a as can be seen from (2.17). From 
Proposition l(ii) we have the derivative 

00 1 lnfrj -I- r/^l 

7 J(a) = C'(2, a) + C(2, a) = C'(2, a) + ^(o) = £ , I m ; » ( 2 - 33 ) 

where the term £(2, a) = V^* 3 ) i s the trigamma function. The function 7 J(a) — > 00 
as a —>■ + and — > through negative values as a — > 00. The sole zero of 7l (a) occurs 
at oj, where 7i(a^) > 0. Therefore, 71 has the single global maximum as claimed. 
Remark. The approximate value 71 (a*) ~ 0.0379557. 



dt 

N+1/2 (t + a) s+1 
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Proposition 6. As described below, we have for Re s < 1 



( {j \s,a)da = 0. (2.34) 



o 



We differentiate (2.8) j times with respect to s and apply this equation, putting 
z — 1 — s, giving the Proposition. 

That (2.34) holds can be seen by first noting that by the functional equation of 
the Hurwitz zeta function we have for Re s < 1 and < a < 1 

C(s, a) = 2 s vr s - 1 r(l - s) £ sin (Wna + ^) n s ~\ (2.35) 

n=l ^ ^ / 

By using the product rule, various forms of £^'(s, a) follow, and these give (2.34). 

Remarks. Equation (2.34) could also be found on the basis of various integral 
representations for ( holding for Re s < 1. 

Proposition 6 gives a generalization of Proposition 4 of [7] when j > 0. 

As a byproduct of our proof of Proposition 6 we have 
Corollary 4. We have 

E 2 ^f^ = (-l) n [C (n) (0,a)+n!]. (2.36) 

k=0 

This is an extension of the "beautiful sum" at a = 1 attributed to O. Marichev [33J. 

As a special case, we have 
Corollary 5. We have 

E = -C'(0, a) - 1 = \ ln(27r) - In T(a) - 1. (2.37) 
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Indeed, Corollary 4 is consistent with the relation (cf. [7], p. 610) Rj(a 
.(_l)j£0-)(0, a ), where 

In-* (n + a) In- 7 a 



— =: — = -7,- In 7 a - > 



n + a 



J>0. 



(2.38) 



By integrating this equation we have 

oo 

Rj(a) - Rj(l) =j(l- a)7j-i - ln J a - ^ 



n=l I 



m J (n + a) —j (a — 1 ) 



r? 



(2.39) 



We record an expression for (— l) J ' +1 i2j(a) in the following. 
Lemma 2. We have for integers j > 1 



C«\s, a) = (-l)V- f ({) (j - k)\ {a ^*_ k+1 + { -^a- In' a 



fc=0 



giving 



00 Pi(x) 



o (x + a) 



s+l 



ln J L (x + a) dx- {-I) 3 s 



00 PAx) 



o (x + a) 



s+l 



W(x + a) dx, (2.40) 



C W (0,a)=a^f^(j-fc)!(-l) fe+1 ln fc a+^ln^+(-iyj h ''' ' 



fc=0 



a X 



Pi(x— a)dx. 
(2.41) 



Lemma 2 follows from a direct calculation using the integral representation [36J (2.3) 
valid for Re s > — 1, 



C(s,a) = — 



a s a 1 ' 



<-°° PAx) 
s I - : — —dx. 



2 s — 1 Jo (x + a 



\s+l 



(2.42) 



Equation (2.40) may also be proved by induction. The a = 1 reduction of (2.42) is 
well known [32J (p. 14). 
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Proposition 7. We have from [36] (pp. 153-154) 

f°° ln n_1 x 
C n (a) = / PAx — a) — On — \nx)dx. 

Ji rrl 



x z 



(2.43) 



We put e n = exp(n/2). We note that the generic function \n n x/x 2 is nonnegative 
and monotonically increasing for 1 < x < e n and nonnegative and monotonically 
decreasing for e n < x < oo. We have that Pi(x — a) is bounded and integrable and 
we split the integrals in (2.43) at e n -\ and e n . By the second mean value theorem for 
integrals (e.g., [16J, pp. 1097-98) we have 



en ln n x ln n e n f en 

— — P\(x — a)dx = — - — / P\(x — a)dx, 
i x 2 e z n J v 

for some 77 with 1 < 77 < e n . Therefore, we obtain 

ln n x 



(2.44) 



1 ar 



-Pi(x — a)dx 



n 



1 n n 

PAx — a)dx < . 

2 n e n J v n ; ~ Q2 n e n 



(2.45) 



Here we have used the standard Fourier series for P\(x) (p. 151) or [TJ (p. 805), 

sin(2n7rx) 



n=l 



TITY 



giving 



so that 



r 1 ~ cos[2mr(x-a)] 

/ P 1 {x-a)dx = -\^ — 

Jb 2 ~ x n z 7i 2 



/ P±(x — a)dx 

Jb 



5 i« 2) = 5- 



(2.46) 



(2.47) 



(2.48) 



Similarly, we have 



00 In™ x In™ e ft 

— — P\{x — a)dx = — y~ I Pi( x ~ o)dx, 



(2.49) 
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for some £ with e n < £ < oo. This gives 

r°° ln n x . , n n ft , 1 n n . 

\ —z-Px(x - a)dx = / PAx-a)dx< . 2.50 

Je n x 2 K ' 2 n e n Je n V ; ~ 6 2 n e n K ' 

Combining the 4 integral contributions of (2.43) yields the Proposition, as 



\Cn{a)\<\ 



(n — I) 11 1 n n n 
+ 



en 

< . 2.51 

- 2 n e n K ' 



2n-l e n-l ' 2 n e 

Remarks. Zhang and Williams previously found the better bound 



However, our presentation shows possibility for improvement. If, for instance, a 
better bound is found for the integrals J 6 C P\{x — a)dx, then we may expect a tighter 
estimation. 

Let us note also that the use of (2.40) for C (pk + 1, a) together with the expres- 
sions of Proposition 2 provides many other opportunities for the estimation of |7fc(a)| 
and |Cfc(a)|. 

Proposition 8. We apply the formula of Klusch [22] (2.5), 

L(x, s, a + = £(-l) fc UiL(x, s + k, a)^. (2.53) 

This formula may be obtained by Taylor expansion, expansion of an integral repre- 
sentation, or by binomial expansion in (1.15). We take n derivatives of (2.53) using 
the product rule. We evaluate at s — 1 using Lemma 1 and the defining relation 
(1.17) of l n . Separating the k = term yields the Proposition and then Corollary 3. 
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Remark. As very special cases, we have 

1 + a 



and 



2 -,a =- ln 2 2 



1 + a 



2 

1 + a 



+ 7i 



1 + a 



71 1.2 



+72 ( ^ ) - 72 



+ 21n2 

a + 1 



7i I 2 I - " 



a + 1 



(2.54) 



As another example, from 



£ 

n=0 



n + a) 



4- C 



, a + 2 
((s,——)+i 



CI 8, 



a + 1 



c u 



a + 3 



(2.55) 



by taking the real and imaginary parts of derivatives evaluated at s — 1, we obtain 
expressions for j k (a/A) — jk[( a + 2)/4] and 7fe[(a + l)/4] — 7fc[(a + 3)/4]. Evidently 
for x a rational number such relations always exist. 

Proposition 9. From the classical theta function-based representation of the Rie- 
mann zeta function [6] (2) we have 



r(*/2)C(*) = + £ n-T (|,7rn 2 ) + tt-V* g n -i r ( i 7 2 - 1 ,™ 2 ) . (2.56) 

We expand both sides of this equation in powers of s — 1, and equate the coeffi- 
cients of (s — 1)° and (s — l) 1 on both sides to obtain parts (i) and (ii), respec- 
tively. The coefficient of (s — l) n_1 for n > of the term 7r s / 2 /[s(s — 1)] is given by 
7r 1 / 2 (— l) n Z)j =0 (— 1) J In 7 7r/2 J , and the simple polar term from the left side of (2.56) is 
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cancelled by the n = term. We use the special function relations T(0, z) = — Ei(— z) 
and r(l/2, z) = \/7r[l — erf(-^/i)]. The incomplete Gamma function [16J is given for 
Re x > by 



r(a, x) 



oo 2x a e~ x f°° t 1 ~ 2a e 
e-H a - l dt = — r / — dt. 



(2.57) 



T(l -a) Jo t 2 + x 

where the latter form holds for Re a < 1. It is convenient in finding derivatives of 



this function to use the relation 



T(a, x) = T(a) lFAa; a + 1; —x), 

a 



a £ N, 



(2.58) 



where iF\ is the confluent hypergeometric function [HUB]- It follows that 



— iF%(a; a + 1; —x) = — [iFi(a; a + 1; —a;) — 2-^2 (ofj a; a + 1, a + 1; —a;)] , (2.59) 
aa a 

and 

G? X a 

—T(a, x) = r(a)i/;(a)-\ — - [—a lnx i-F 1 (a; a + 1; — x) + 2-^2(0:, a; a + 1, a + 1; — x)] . 

(2.60) 

When a — > in (2.60), the singular terms 1/a 2 cancel. In fact, we have the expansions 

(2.61) 



1 7 2 7T 2 

r'(a) = T(a)V(a) = ^ + 4r + — + 0(a), 



a 2 2 12 



and 



x (a J j 
a 2 (a + l)j 



-a In x + 



a 



^— x) J (1 — jinx) 



(-x) j + 0(a), a^O. 



(2.62) 
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Summing the latter relation on j as a — > then gives 
00 (I — i \nx) 

E -3/ n H J = ~ x 3^3(1, 1, 1; 2, 2, 2; -x) + [7 + T(0, z) + lure] lnx. (2.63) 

Pi ru - !) ! 

We apply these relations at a = and 1/2 and the Proposition follows. 

Remarks. By making a change of variable in the theta function-based represen- 
tation [6] (2) we may include a parameter in the series representation (2.56), and so 
in Proposition 9 too. Beyond this, we may include a parameter b > and a set of 
polynomials Pj{s) with zeros lying only on the critical line in representing the Rie- 
mann zeta and xi functions [12] (Proposition 3). Thereby, we obtain a generalization 
of Proposition 9. 

In equation (1.20), the sum terms provide a small correction ~ 0.0230957 to 
produce the value 7/2 ~ 0.288607. Although there are infinite sum corrections in 
Proposition 9, the expressions such as (1.20) and (1.21) may have some attraction for 
computation. This is owing to the very fast decrease of the summands with n. For 
(1.20), using known asymptotic forms, the summand terms have exponential decrease 
~ g-" 2 ^ _|_ q _ Similarly for (1.21), the summand has exponential decrease 
in n. Therefore, high order approximations for the constants may be obtained with 
relatively few terms. 

We note an integral representation for the term — Y^=i Ei(— 7m 2 ) in Proposition 
9 in the following. 

Lemma 3. Put the function 9^{y) = 1 + 2X^i2/™ 2 - Then we have 

dx 



-£Ei(-™ 2 ) 1/ e7r /(,-i) 

n=l ^ 



X — 1 
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1 roc An, 

Proof. Let L° be the Laguerre polynomial of degree n and parameter a. We use 
the relation (e.g., [IS], p. 1038) 

at a = 0. We then have 



r(a,x) = x a e- a; 5:-^, «>-!, (2.65) 



- E Ei(-7rn 2 ) = E e ~™* E I/ '' (7r " 



Z- _i_ 1 

n=l n=l fc=0 ^ T X 



OO OO r.\ oo 

E e ~™ 2 E / x fc dxE^(™ 2 ) 

n=l fc=0 fc=0 



,2 



oo „i xttti 2 /(x— 1) 

= - V e"™ 2 / da;, (2.66) 

where we used the generating function of the Laguerre polynomials (e.g., [16], p. 
1038). The interchange of summation and integration is justified by the absolute 
convergence of the integral. Using the definition of #3 completes the Lemma. 

Discussion: Generalization of Proposition 9 

We may generalize Proposition 9 to the context of generalized Stieltjes constants 
coming from the evaluation of Dirichlet L series derivatives at s = 1. For this we 
require a number of definitions. We use [20] (Ch. 1, Section 4.2). 

Let x be a primitive character modulo k. We need two theta functions, 

00 

2 , 



0( T ,X) = E X(n)e 



-7TTTL j k 
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for x an even character, 

oo 
n=—oo 

for x an odd character, and the Gauss sum 

3W = Ex(j)e 2njA 
i=i 

We let 5 = or 1 depending upon whether x(— 1) = 1 or — 1, respectively. We put 

as,x) = (*k-r {s+1)/2 r (^) L( 8 ,x). 

Then by [20] (p- 15) we find for X an even character, 

2£(«, X) = jH r s/2 -^(s, X )dr + ^ r-(' +1 >/ 2 0(r, x)rfr 



n=— oo 



2" 



+ _^L 7r (-D/2 A; ( 1 -)/ a £ xW^r 



7T- 



n^oo \ 2 ' k 

and for x an odd character, 



7r - (s+ l)/2 A ,( s+ l)/2 g nx ( n ) n -( S +l) r /i±l + 1)7r ^ 



n=— oo 



Then one may multiply the expressions for £(s, x) by (7rfc _1 )^ +1 ^ 2 and proceed as in 
the proof of Proposition 9 in developing both sides in powers of s — 1. 
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In the case of the Hurwitz zeta function we employ the theta function 

6( T ,a) = Y j e-^ n+a) \ 

with its functional equation 

9(l/r,a) = ^e-™ 2/T 6{T,-ia/T). 

Proposition 10. We make use of expressions in the fine paper of Fine [T5j. In 
particular, see his expressions for the function H(a, s) (pp. 362-363). We have 

("1 c \ oo poo 
[C(l -8,a) + C(l - s, 1 - a)} = 2 y cos(27rna) / e'^H^dt 
2 ' „=i 7 ° 



2n- s/2 T (-) y cos ( 27ma ) 



2/ n* 

= vr- s/2 r [Li s (e 2 ™) + Li s (e" 2 ™)]. (2.67) 

Initially, the left side of this equation is valid for Re s < 0. With analytic continuation, 
the function is H(a, s) + 2/s is entire in s. Therefore, we may expand both sides of 
the representation (2.64) about s = 0. As it must, the singular term —2/s effectively 
cancels from both sides of (2.67). Part (i) of the Proposition results from the constant 
term s° on both sides of (2.67), and using the reflection formula for the digamma 
function — a) = ip(a) + 7rcot7ra. Similarly, part (ii) results from the term s 1 on 
both sides of (2.60). 

Remarks. Proposition 10 and (2.67) properly reduce as they should for a = 1/2. 
In this case, we have the alternating zeta function J2^=i = (2 1_s — 1)C( S )- Then 
both sides of part (i) of Proposition 10 yield 7 — In it + 2 In 2. 
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Proposition 11. By using the definition (1-24) we have 

1 rv 

Pn+l(y) = —, i X{1- X) n -i dx 

n\ Jo 

n-l 



1 " _1 ry 
— Y(-l) n+k - l s(n-l,k) x(l-x) k dx 



k=0 

- ^g-^-'^gfl, [(% + » + !)(!- V)\V - 1) + 1], (2-68) 

whence part (i) follows. Of the many ways to perform the elementary integral here, 
it may be evaluated as a special case of the incomplete Beta function (e.g., [16], p. 
950) 



1 rv (— l) n+1 n fy 

Pn+M = 7 (-z)n dx = 1 \ £ s{n, k) / x k dx. (2.70) 

n! Jo n\ r-L Jo 



B x (p,q) = [ X t p - 1 (l-t) g - 1 dt = — 2 F l (p,l-q;p+l;x), Re p > 0, Re g > 0. 
jo p 

(2.69) 

For part (ii) we again use a generating function relation for s(n, k), writing 

1 fy, , . f-l) n+1 " 

fc=0 

The second line of (1.26) simply follows from the values s(n, 0) = 5 n o and s(n, n) = 1. 
The identity of part (hi) follows simply from putting y = 1 in parts (i) and (ii). 
For part (iv), we have the well known expression 

7fc = lim ( V — - - J— ln fc+1 n) , (2.71) 

so that lim^oo rffi = 0. We form 

«jf=f: f— - f +1 — 



m=l \ m Jm X 
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n rl 

E 



m=l 



\n k m \n k (x + m) 



rn 



x + m 



dx. 



(2.72) 



Therefore, we have 



oo r \ 



E 

m=n+l 



r hr fc m ln fc (:r + m) 



m 



x + m 



dx. 



(2.73) 



We now use pDj (5.7) 



\n n y \n n (x + y) ^ x k 



x + y 



k=i y 



i n—l 

(-i) fe ^i/ + 7jE 



nl 



s{k + l,j + 2) In^^y 
(2.74) ' 



Performing the integration of (2.73) gives part (iv). 



Discussion related to Proposition 11 

By defining the function fk(x,m) = (x + m)\n k m — m\n k (x + m), it is possible 



to further decompose the remainders of (1.29) by writing 



oo „i 

E / fhix,m) 

=n+l J ° 



m(x + m) m{m + 1) 



dx+ V 



m=n+l 



fk(x,m) 
m(m + 1) 



(2.75) 



The representation (2.74) can then be used three times in this equation. This process 
of adding and subtracting terms can be continued, building an integral term with 
denominator m(m + 1) • • • (m + j)(m + x). For k = 0, there is drastic reduction to 
the original construction of Ser [29] . 

Proposition 11 (iv) can be easily extended to expressions for r^ k \a) = 7fc(a) — 

As pointed out by Ser, it is possible to write many summations and generating 
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function relations with the polynomials P n . As simple examples, we have 

E K + i(y)z n = - E (-1)" f y ) z B = 1 - (1 - *)», (2.76) 



n=l n=l 

and 



co r>i 



E = - E — (! ] = + 1) + 7 = H y . (2.77) 



n=i n n \ n 



With the form (1.27) we easily verify a generating function: 

En y n-i _ (-1)" +1 -A s(n,k) ^ n _ x 

n=l n=l n! fc=l K + 1 

1 ' ] < ^ ln fc (l-^ 



(*+!)« 

= - + wi — 2 - 78 

z ln(l — 

One may also write a great many hypergeometric series of the form J2™=i — ^rj — 2,1 > 
that we omit. 

In contrast to an expression of Ser [29], we have 

oo 111 
^ 2 1-e 2 2 

11 , /2 

= cotn - 

z 2 V2 

CO n 

where the latter series may be found from [IB] (p. 35). This relation serves to connect 
the numbers p n +i with the values C(2&). Equivalently, using an integral representation 
for C(2/c), we have 

EPn + i( 1 -e 2 r 1 = 2/ sinhM^— -. (2.80) 

n=2 J0 e 1 
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As well, we have many extensions including 

OO I tz 1 

- cT-^-V = |n | 1 + eJ(1 _„_ e _„i + - 5. (2.81) 

We may inquire as to the asymptotic form of p n +i as n — ► oo. Using only the 
leading term of the result of [31], we have 



1 " (-l) fc In^n 1 
^^ng l ( Hl)(^l)! ^' n ^°°- (2 ' 82) 

There are corrections to this result which overall decrease it. Therefore, one may 
wonder if it can be proved that the expression on the right serves as an upper bound 
for all > 2. 

Motivated by pQ (p. 824), we conjectured that 

1 

Pn+i 7] 1 \o> n -> °°> ( 2 - 83 ) 

is an improved asymptotic form and upper bound. This asymptotic form has been 
verified by Knessl [24]. In fact, he has obtained a full asymptotic series for p n +i in 
the form 

1 / 00 A \ 

i^i~-tt- 1 + £rr- ' (2 - 84) 

n in n \ J=1 hr n J 
where Aj are constants that may be explicitly determined from certain logarithmic- 
exponential integrals. We have A\ = —27 and A 2 = 37 2 — tt 2 /2. Knessl develops 
(2.84) from the exact representation 

P** 1 = / 7TZ wi 2 T iT » n ^ L ( 2 - 85 ) 
Jo (1 + u) n (In u + n z ) 
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From (2.85) we have the following 
Corollary 6. We have the integral representation 



7= f e*Mi±£_2 d2 . 



Proof. We have pj] (p. 943) 



7 = -^(1) = + T^—) du = I' ( r~7T~ 7 + - J d« 

jo Vmu 1 — u/ Jo \m{l — v) v) 



= ffpn + i^ = E— > ( 2 -87) 
where we used (2.78). Now we use Knessl's representation (2.85), giving 

™lr°°l du 
7 = L 



^nJo (1 + u) n (In 2 u + 7i 2 ) 



°° 1 / ?i \ G?1t 
m 



Vl + ^7 (In 2 U + 7T 2 ) 

/OO q z 
_ 7^T^t [ " " 1,1(1 + e ' )]<iz 

ln(l + e" 2 )^. (2.87) 



-oo \Z 2 + 7T 2 ) 

oo e z 



-oo (z 2 + 7T 2 ) 

Summary and very brief discussion 

Our methods apply to a wide range of functions of interest to special function 
theory and analytic number theory, including but not limited to, the Lerch zeta 
function and Dirichlet L functions. Two dimensional extensions would be to Epstein 
and double zeta functions. Our results include an addition formula for the Stieltjes 
coefficients, as well as expressions for their derivatives. The series representations for 
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7fc(o) have very rapidly convergent forms, making them applicable for multiprecision 
computation. Byproducts of our results include some summation relations for the 
Stieltjes coefficients. We have given a means for estimating the magnitude of these 
coefficients, although this remains an outstanding problem. 
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